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FOREWORD 


Mr.  Thompson  died  suddenly  in  1970  with  this  report  neadv  completed.  Shortly  thereafter,  SFC 
Scherer  separated  from  the  Army.  The  manuscript  was  gradually  forgotten  until  rediscovered  in  a 
cleanup  in  1988.  The  subject  of  acoustic  location  is  of  ongoing  interest  to  the  Army  (sec  Field 
Artillery  Sound  Ranging,  FM6-122,  December  1983). 

Acoustics  research  at  Aberdeen  Proving  Ground  no  longer  includes  location  and  ranging  projects 
to  the  editor’s  knowledge.  The  topics  arc  vehicle  acoustic  signatures,  muzzle  blast,  noise  relation  to 
troop  safety,  and  especially  community  (upper  Chesapeake  Bay)  noise  exposure.  The  forward 
detectors  the  author  speculated  on  for  accurate  location  of  artillery  are  not  in  use;  but  unattended 
sound  level  recorders  are  in  routine  use  by  Range  Control,  Combat  Systems  Test  Activity  for  noise 
monitoring. 

An  acoustic  focusing  model  is  run  at  BRL  for  predicting  noise  levels  at  a  location  where  damage 
claims  were  made  against  the  Army. 

There  is  independent  reason  to  suppose  that  accurate  acoustic  location  of  artillery  may  be 
possible.  Naturalists'  studies  of  nocturnal  predators  (bam  owls)  show  that  the  predator  depends  on 
listening  to  locate  and  strike  its  prey. 
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1.  INTRODUCTION 


The  Ballistic  Research  Laboratory  (BRL)  was  requested  to  investigate  the  problem  of  delecting 
and  locating  certain  foreign  tactical  artillery  rockets  by  acoustic  techniques.  The  rockets  of  interest 
include  107-mm,  122-mm,  and  140-mm  projectiles.  The  Soviet  140-mm  artillery  rocket  (Model 
M-14-OF)  is  shown  in  Figure  la.  The  122-mm,  single  tube,  "6-ft"  rocket  is  shown  in  Figure  lb. 

The  Chinese  107-mm,  Type  63  rocket  is  shown  in  Figure  lc. 

Research  at  BRL  has  produced  models  for  predicting  the  field  propagated  by  a  variety  of 
stationary  and  moving  sources  of  acoustic  energy.  This  report  gives  the  model  prediction  of  noise 
around  the  Chesapeake  Bay  from  the  largest  rocket.  This  example  should  provide  some  indication  of 
from  where  and  how  well  die  launch  position  of  these  rockets  can  be  located  by  acoustic  techniques. 

2.  THEORY 

2.1  Prediction  Models.  Many  theoretical  and  empirical  assumptions  were  required  to  make  a 
mathematical  model  which  could  be  solved  for  particular  applications.  A  pan  of  die  BRL  mode!  is 
based  on  the  acoustic  wave  equation  (die  ray  physics,  high  frequency  approximation  in  punicular). 
This  is  complemented  by  scmicmpirical  rclauonships  in  areas  where  the  wave  equation  or  its 
approximation  are  difficult  to  apply.  The  appropriate  equations  of  ray  physics  arc  used  for  calculating 
die  position,  density,  and  time  sequence  of  acoustic  rays  in  the  layered  atmosphere,  while 
scmicmpirical  relationships  arc  used  for  intensity  predictions  that  depend  on  both  die  source  and 
atmosphcnc  absorption  characteristics.  In  general,  we  employ  more  than  one  mediod  for  intensity 
spectrum  predictions  to  allow  comparison  between  results  obtained  using  different  assumptions. 

A  brief  descripuon  of  semicmpincal  relationships  used  with  the  acoustic  ray  physics  approximation 
is  given  in  the  Appendix.  Further  details  on  the  theory  and  applicauon  can  be  found  in  Thompson 
(1968);  Lassiter  and  Heitkottcr  (1954);  Mayes,  Lanford,  and  Hubbard  (1959);  and  Humphrey  (1957). 
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2.2  Detectability  of  Rockets.  Since  the  rocket  motor  noise  and  possibly  the  motor  initiation  or 
cutoff  sounds  arc  the  most  consistently  available  signals  for  launch  point  location  within  a  range  of 
10  km,  these  signals  arc  predicted  down  to  a  magnitude  of  50  dB.  Below  this  level,  ambient  noise 
makes  detection  difficult. 

The  rocket  motor  initiation  signal,  which  contains  frequencies  below  50  Hz,  should  be  more 
detectable  a»  the  greater  distances  than  the  rocket  motor  noise.  In  addition  to  experiencing  less 
attenuation,  these  lower  frequencies  arc  not  as  readily  refracted  away  from  the  ground  by  a  negative 
sonic  gradient. 

The  ballistic  wave  signal,  caused  by  the  shock  in  front  of  a  supersonic  missile,  is  detectable  the 
farthest,  depending  on  circumstances  of  missile  trajectory  and  meteorological  conditions. 

These  three  rocket  signal  categories  (initiation,  motor,  and  ballistic)  can  be  identified  through  the 
use  of  a  spectrogram  showing  the  amplitude  spcctnim  vs.  time. 

2.3  Location  Techniques.  In  order  to  locate  the  launch  position  of  artillery  rockets,  we  can  use 
one  or  both  of  the  following  techniques: 

(1)  Measure  the  arrival  lime  of  the  same  part  of  the  signal  at  three  or  more  positions.  With  this 
information,  techniques  (Dobrin  1976  and  Small  1953)  arc  available  which  can  be  used  for 
locating  the  origin  of  the  signals. 

(2)  Determine  two  or  more  azimuths  along  which  the  signal  is  propagating.  With  this 
information,  the  origin  can  be  located. 

For  applying  the  first  technique,  a  fast  initial  rise  time  of  the  acoustic  signal  is  desirable  to 
facilitate  measurement  of  the  same  part  of  the  signal  at  all  stations.  A  detectable  rocket  motor 
initiation  signal  would  permit  use  of  the  arrival  time  technique,  provided  that  it  could  be  distinguished 
from  the  (sometimes  present)  ballistic  wave  signal. 
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The  second  technique  is  described  in  sound  ranging  literature  (Fox  1967a;  Fox  1967b;  Bellucci 
1966;  Barichivich  1966;  and  Nordquist  1967).  For  applying  the  second  technique,  each  detection 
station  must  have  a  small  array  of  sensors  arranged  so  that  their  outputs  can  be  correlated  to  determine 
the  direction  of  propagation  of  a  signal  across  the  array. 

Ballistic  wave  signals  could  be  detected  by  certain  station  arrays,  and  continuous  rocket  motor 
noise  (or  motor  initiation  signals)  could  be  detected  by  other  arrays.  A  small  sensor  array  has  been 
suggested  to  determine  azimuths  of  acoustic  propagation,  since  it  has  been  proven  successful.  Other 
methods  of  determining  azimuth  should  also  be  considered,  such  as  the  use  of  vector  particle  velocities 
and  pressures  measured  together.  Vector  particle  velocity  measurements  have  recently  become  more 
feasible. 

From  the  above  discussion,  we  can  see  that  a  practical  technique  for  location  depends  on  the 
nature  and  detectability  of  the  signal.  Whether  either  technique  can  be  employed  must  be  determined 
from  feasibility  tests  designed  from  prediction  models  such  as  those  of  Thompson  (1968). 

The  accuracy  of  launch  point  location  for  the  First  (arrival  time)  technique  would  depend  on  (a) 
the  precision  of  the  amval  time  of  the  phase-correlated  signal  measured  at  the  stations,  (b)  how  well 
die  sonic  velocities  between  the  source  and  the  different  stations  arc  known  (this  factor  becomes 
significant  at  detector  distances  greater  than  2  km),  and  (c)  the  distances  between  the  stations  and  the 
source.  The  accuracy  obtained  using  the  second  (azimuth)  technique  would  depend  mostly  on  the 
precision  of  the  measured  angles  and  the  source -detector  distances. 
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3.  PROCEDURE 


Predictions  in  this  report  are  presented  for  a  solid  fuel  rocket  similar  to  the  foreign  140-mm 
artillery  rocket  mortar  (Model  M-14-OF;  Figure  la).  Both  spectral  and  peak  intensities  were  estimated 
at  distances  considered  reasonable  for  Held  location  of  the  launch  point. 

A  static  test  was  assumed  in  a  homogeneous  atmosphere  for  the  initial  predictions.  The  results 
obtained  under  this  assumption  provided  some  basic  understanding  of  the  propagation  of  rocket  motor 
initiation  signal  and  the  rocket  motor  noise  without  the  complications  of  weather  and  flight  path. 

Since  some  acoustic  measurements  were  taken  by  the  Air  Force  (Wathen-Dunn,  unpublished)  at  about 
100  m  from  this  foreign  rocket  Fired  near  Edgewood  Arsenal,  MD,  this  measured  spectrum  was 
projected  to  greater  distances  using  the  BRL  model,  then  compared  with  the  predictions  for  the  U.S. 
solid  fuel  rocket  having  similar  characteristics. 

For  dynamic  test  predictions,  a  set  of  meteorological  parameters  (with  temperature  inversion)  is 
specified,  since  winds  and  temperature  should  be  considered  in  an  accurate  source-locating  technique. 
For  discussion  purposes,  the  meteorological  condition  existing  on  26  October  1964  at  Aberdeen 
Proving  Ground,  MD,  was  selected.  This  condition  corresponds  to  that  existing  for  Case  2  in 
BRL-MR-1930. 


Predictions  for  dynamic  tests  are  given  as  maximum  surface  intensity  levels  for  both  the  rocket 
motor  noise  and  ballistic  energy  signals.  Although  the  directional  ballistic  energy  signal  can  be 
detected  at  greater  distances  along  the  surface  under  selected  conditions,  the  rocket  motor  initiation 
and  rocket  motor  burning  signals  can  be  detected  under  most  prevailing  conditions  and  in  most 
directions. 

With  the  above  prediction  estimates  under  the  specified  conditions,  some  results  and  conclusions 
arc  given  about  detection  and  location  of  this  type  of  rocket.  In  discussing  the  results,  both  the  arrival 
lime  and  azimuth  techniques  for  launch  point  location  are  considered. 
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4.  RESULTS 


Figure  2  shows  intensity  spectra  (decibels  of  sound  vs.  frequency)  predicted  at  three  distances  on 
the  surface  in  a  homogeneous  atmosphere  from  a  statically  Fired  U.S.  solid  fuel  rocket  similar  in  size 
and  thrust  to  the  foreign  140-mm  rocket. 

Figure  3  shows  the  predicted  peak  intensity  vs.  distance  curve  (solid  line)  for  the  U.S.  rocket. 

Also  shown  in  Figure  3  (dashed  line)  is  the  intensity  vs.  distance  curve  produced  by  spherical 
divergence  of  the  wave  front  alone,  i.e.,  no  atmosphere  frequency  absorption. 

Figure  4  shows  a  comparison  of  the  10-km  spectrum  in  Figure  2  with  a  10-km  spectrum  deduced 
from  Air  Force  measurements  made  100  m  from  the  foreign  artillery  rocket  when  it  was  launched  at 
Edgewood  Arsenal  (Wathen-Dunn,  unpublished). 

Figure  5  shows  the  weather  conditions  chosen  for  the  dynamic  noise  predictions. 

Figure  6  shows  the  maximum  surface  signal  intensity  contours  predicted  for  the  foreign  solid  fuel 
rocket  motor  if  it  were  launched  easterly  over  the  Chesapeake  Bay  under  the  weather  conditions  shown 
in  Figure  5. 

Figure  7  shows  the  maximum  ballistic  wave  signal  intensities  that  occasionally  can  reach  the 
surface  if  the  140-mm  rocket  is  fired  at  an  azimuth  of  89°,  at  a  low  angle  (3°)  under  the  weather 
conditions  in  Figure  5.  It  should  be  emphasized  that  Figure  7  shows  only  ballistic  wave  signals 
coming  from  the  rocket.  All  the  other  figures  show  only  the  rocket  motor  initiation  and  burning 
signals. 
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Intensity  vs.  Distance  Curve  for  Solid  Fuel  Rocket  (3,500-lb  Thrust 
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TEMPERATURE 


F'li'urc  7  Predicted  Surface  dH  Contours  for  the  Ballistic  Wave  Energy  From  the  140-mm 
Rocket  Artilicrv  Fired  at  a  Low  Ancle  Dunne  the  Selected  Weather  Conditions 


5.  DISCUSSION  Oh  RESULTS 


5.1  Detection.  F-'rom  the  frequency  spectra  in  Figure  2,  we  can  see  that  detection  of  artillery 
rockets  having  about  3.500  lb  of  thrust  is  possible  for  distances  greater  titan  10  km  under  low  ambient 
noise  levels.  Higher  background  noise  levels  from  moderate  winds  would  severely  reduce  this  range 
of  detectability,  unless  wind  screens,  filters,  or  other  noise  reduction  techniques  were  used. 

Other  factors  influencing  signal  detectability  that  were  not  taken  into  account  in  Figure  2  include 
weather  conditions,  height  of  rocket  motor  cutoff,  propellant  temperature,  and  acoustic  directivity  of 

the  source. 

In  Figure  4,  the  l()-km  spectra  of  the  U  S.  and  foreign  motor  look  somewhat  similar.  It  is 
at  these  greater  distances,  of  course,  that  identification  and  location  of  source  arc  most  needed. 

The  foreign  spectrum  is  a  distance  projection  from  the  BRL  model  and  is  based  on  only  one  set  of 
measurements  made  very  close  ( 100  mi  to  the  trajectory.  The  projection  uses  divergence  and 
absorption  laws  formulated  like  Equations  A-3,  A-4,  and  A-5  in  the  Appendix. 

Other  factors  besides  divergence  and  absorption  must  be  considered  before  w-c  can  effectively 
project  these  very  close-in  measurements  to  the  distant  environment.  Some  examples  of  these  factors 
will  now  briefly  be  discussed. 

The  measurements  (Wathen-Dunn,  unpublished)  that  we  projected  (Figures  6  and  7)  were  made 
very  close  to  a  moving  source  that  generated  Doppler  effects  w'hich  changed  with  time,  direction,  and 
distance.  A  gage  was  located  where  leaking  (indirect)  ballistic  wave  signals  could  be  recorded  along 
with  rocket  motor  initiation  and  rocket  motor  burning  signals.  The  more  directional  ballistic  wave 
signal  would  have  considerably  more  variation  with  location  and  trajectory  than  the  less  directional 
rocket  motor  and  staging  signals.  The  response  of  the  measuring  sensor  dropped  off  rapidly  below 
30  Hz,  where  initiation  signals  may  have  substantial  energy.  In  spite  of  these  problems,  it  is  fortunate 
that  the  Air  Force  measurements  were  made. 

The  dashed  line  in  Figure  4  shows  that  the  most  suitable  frequencies  for  detection  at  10  km 
arc  40  and  2(X)  Hz.  If  the  Air  Force  instruments  could  register  below  40  Hz,  it  is  expected  that  the 
40-11/  peak  shown  in  the  dashed  line  in  Figure  4  would  be  shifted  to  a  higher  peak  near  20  Hz. 
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This  expected  low  frequency  peak  is  attributed  to  the  rocket  motor  initiation  pulse,  while  the 
200-Hz  peak  is  attributed  to  rocket  motor  noise. 

It  is  suggested  that  signal  recordings  to  a  range  of  10  km  be  obtained  for  revised  predictions. 
These  measurements  should  have  flat  frequency  response  from  10  to  500  Hz. 

Spectrograms  similar  to  those  presented  in  the  Air  Force  report  (Wathen-Dunn,  unpublished), 
but  made  from  measurements  at  greater  distances  (1-10  km)  and  including  a  lower  frequency  range 
(to  5  Hz),  are  suggested  to  explore  rocket  identification  from  acoustic  signals.  Preferably,  a  spectral 
contour  plotter  would  be  used  to  present  the  signals  obtained  at  this  1-km  to  10-km  range.  These 
contour  plotters  arc  commercially  available.  Once  the  signature  has  been  measured  and  its  energy 
spectrum  analyzed,  a  readily  available  artificial  source  (Galloway,  Watters,  and  Baruch  1965)  could 
be  found  to  test  various  location  techniques  under  different  environmental  conditions. 

The  peak  intensity  levels  from  the  burning  motor  were  estimated  at  intermediate  distances  with  a 
peak  intensity  vs.  distance  curve  similar  to  that  shown  in  Figure  3.  The  difference  between  the  dashed 
curve  (showing  only  spherical  divergence)  and  the  lower  solid  curve  (showing  the  predictions)  results 
from  acoustic  absorption,  which  increases  with  frequency.  Since  theoretical  absorption  laws  do  not 
agree  with  absorption  as  measured  from  rocket  motors,  semiempirical  laws  for  absorption  were 
selected  in  preference  to  the  purely  theoretical  estimates  (see  Equation  A^t  in  the  Appendix). 

In  order  to  obtain  realistic  estimates  of  the  possible  acoustic  signals  for  a  dynamic  rocket  flight 
under  typical  weather  conditions,  we  would  have  to  present  data  for  several  weather  and  flight 
conditions.  Since  this  report  is  not  written  for  such  extensive  treatment,  one  weather  condition  was 
selected  for  maximum  intensity  predictions  resulting  from  both  the  rocket  motor  signals  and  the 
ballistic  wave  signals. 

From  Figure  6,  which  shows  the  distribution  of  acoustic  intensity  from  the  rocket  motor  we  can 
sec  'hat  the  60-dB  contour  extends  more  than  10  km  in  the  southerly  direction  and  less  than  10  km 
in  the  northerly  direction.  From  the  weather  profile  in  Figure  5  >'  can  be  deduced  that  this  lack  of 
acoustic  symmetry  in  Figure  6  results  from  a  sonic  inversion  of  the  south  and  the  absence  of  such  an 
inversion  to  the  north.  Although  a  60-dB  signal  is  detectable  over  the  ambient  noise  level  for  most 
'tight  wind  conditions,  more  signal  amplitude  may  be  necessary  for  sufficient  phase  correlation  for  use 
of  the  arrival  time  location  technique. 
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Figure  7  shows  that  ballistic  wave  signals,  much  stronger  than  the  motor  signals,  are  generated  if 
the  rocket  exceeds  the  speed  of  sound.  However,  these  signals  reach  the  surface  only  in  the  direction 
of  launch  (89°  in  Figure  7),  and  only  for  very  low  firing  angles,  and  only  if  sonic  inversions  exist  in 
the  launch  direction.  The  possible  existence  of  a  ballistic  wave  signal  must  be  considered  since  its 
earlier  arrival  may  give  an  erroneous  location  if  used  in  the  arrival  time  technique  with  the  rocket 
motor  initiation  signal  obtained  at  another  station. 

5.2  Location  Accuracy.  Although  accuracy  of  launch  point  location  was  not  considered  in  the 
results  in  this  report,  it  must  be  discussed.  It  is  difficult  to  estimate  the  possible  accuracy  without 
studying  the  nature  of  the  acoustic  signals  measured  at  distances  between  3  and  10  km.  Such  signals 
would  be  used  in  the  initial  attempt  at  rocket  launch  point  location.  From  such  signal  measurements, 
we  must  determine  various  parameters  such  as  phase  correlatability,  frequency  content,  the  motor 
initiation  signals,  rise  times,  and  other  characteristics.  Some  of  these  effects  are  discussed  in  the 
literature  on  sound  ranging  from  guns.  It  is  not  the  purpose  of  this  report  to  discuss  published  sound 
ranging  techniques  that  could  be  partially  applied  to  rocket  artillery.  The  reader  is  referred  to  Fox 
(1967a;  1967b),  Bellucci  (1966),  Barichivich  (1966),  and  Nordquist  (1967). 

From  the  azimuth  and  arrival  time  errors  reasonable  for  these  foreign  rockets,  it  can  be  estimated 
that  a  launch  point  location  error  of  less  than  200  m  is  not  likely  for  such  artillery'  with  measurements 
made  at  distances  between  4  and  10  km. 

A  good  potential  exists  for  improving  the  accuracy  of  acoustic  location  techniques  if  applied 
research  efforts  are  directed  to  this  problem.  With  a  slight  improvement  in  the  state-of-the-art  of  the 
1970s,  acoustic  techniques  should  be  useful  for  identification  and  approximate  location  of  launch  sites. 
This  identification  may  allow  a  more  immediate  notification  of  any  rocket  artillery  attack.  Such 
notification  could  enable  more  people  to  take  cover  and  also  provide  more  time  for  counterattack. 

In  order  to  more  accurately  locate  rockets  that  can  be  used  against  a  community,  an  appropriate 
monitoring  system  must  include  detectors  between  the  launch  point  and  the  community.  This  system 
might  consist  of  a  pattern  of  battery-operated  telemetering  sensors  located  around  the  community. 
These  sensors  would  send  back  any  large  (or  impulsive)  signal  to  a  central  receiver  in  the  community. 
By  timing  the  signal  arrivals  at  several  sensors  close  to  the  source,  the  source  could  be  more 
accurately  located.  Such  a  monitoring  system  is  being  built  for  test  purposes  at  Aberdeen  Proving 


Ground  in  order  to  control  the  problem  of  community  annoyance  from  explosive  sources.  This  gage 
package  is  being  designed  to  telemeter  only  those  signals  exceeding  a  prescribed  intensity  level.  The 
package  will  have  additional  sensors  attached  to  lines  (50-ft  lengths)  that  will  turn  on  telemetry  of  the 
central  sensor  in  time  to  respond  to  a  blast  wave  greater  than  a  selected  intensity  level. 

Another  technique  which  can  be  used  to  improve  accuracy  is  to  fire  a  live  round  back  on  the 
estimated  target  location  and  compare  the  explosion  signal  arrival  times  (or  azimuths)  from  this  round 
with  those  produced  by  the  rocket  launch.  Appropriate  adjustments,  based  on  acousuc  comparisons, 
could  bring  accurate  counterfire  on  the  launch  position. 

Seismic  signatures  should  also  be  considered  for  use  in  locating  a  rocket  launch  position.  It  has 
been  estimated  that  the  seismic  air-induced  signals  could  be  detected  beyond  5  km  from  the  launch 
point 

If  an  accurate  zero  launch  time  could  be  determined  for  use  with  the  acoustic  signal,  much  greater 
location  accuracy  would  be  possible.  Although  electric  or  electromagnetic  signals  generated  at  motor 
ignition  would  provide  such  a  zero  time,  the  evidence  supporting  the  possibility  of  detecting  such 
signals  has  not  been  encouraging. 
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6.  CONCLUSIONS 


There  is  a  paucity  of  measurements  of  acoustic  signals  in  the  1-km  to  10-km  range  from  dynamic 
firings  of  solid  fuel  rocket  artillery  weapons  (Figures  la,  lb,  and  lc).  Without  such  data  to  test 
theoretical  assiur.pLo.o,  ii  is  ditncuit  to  predict  uie  identification  potential  or  the  accuracy  witn  wmen 
such  rocket  launches  could  be  located  with  acoustic  detectors.  The  following  preliminary  conclusioas 
can  be  listed  for  later  experimental  verification. 

(1)  It  is  probable  that  the  rocket  motor  noise  around  300  Hz  can  be  detected  out  to  a  range  of 
10  km  under  low  ambient  noise  conditions. 

(2)  Sound  spectrograms  measured  at  1  to  10  km  from  the  launch  site  should  be  useful 
for  signal  identification. 

(3)  Proof  of  the  existence  of  an  impulsive  motor  initiation  or  cutoff  signal  with  energy  between 
10  and  40  Hz  is  needed.  Having  these  signals  would  improve  the  accuracy  with  which  the 
launch  could  be  located  from  distant  detectors.  These  signals,  at  the  lower  frequencies, 
should  have  greater  amplitudes  beyond  1.5  km  than  the  amplitudes  of  the  rocket  motor  noise. 

(4)  The  ballistic  wave  signal  should  occasionally  be  detectable  on  the  ground  at  distances  greater 
than  10  km  in  the  launch  direction  when  the  artillery  rocket  is  fired  supersonically  at  low 
angles  in  a  direction  where  sonic  inversions  exist  to  bring  this  ballistic  wave  energy  to  the 
surface.  Under  these  circumstances,  ballistic  wave  signal  azimuths  measured  at  the  greater 
distances  may  also  be  used  for  identification  and  location  purposes.  If  an  arrival  time 
technique  is  used,  location  errors  will  result  if  the  earlier  ballistic  wave  signal  arrivals  arc 
mistaken  for  a  rocket  motor  initiation  signal.  The  higher  frequency  content  of  the  ballistic 
wave  signal  could  be  used  to  prevent  this  mistaken  identity. 
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(5)  The  most  logical  methods  for  location  should  use  either  correlated  phase  arrivals  or 
calculated  azimuths  at  two  or  more  stations. 

(6)  Launch  location  with  an  error  between  100  and  500  m  should  be  possible  using  detectors 
at  distances  between  4  and  10  km. 

(7)  For  better  location  accuracy,  one  should  give  more  attention  to  battery-operated  telemetering 
detectors  that  can  be  placed  close  to  the  source. 

(8)  Arrival  times  from  a  launch  could  be  compared  with  those  from  the  explosion  of  a  live  round 
Fired  back  at  the  enemy  launch  position.  With  this  information,  counterbattery  Fire  can  be 
adjusted. 

(9)  Other  signals  (besides  acoustic)  propagating  from  the  launch  site  should  be  further 
investigated.  If  such  signals  exist,  a  better  locating  system  using  several  signal  types  could 
be  developed. 
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APPENDIX: 

SEMIEMPIRICAL  RELATIONSHIPS 
USED  WITH  THE 
RAY  PHYSICS  APPROXIMATION 
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The  ray  physics  approximation  to  the  wave  equation  can  be  very  useful  for  predicting  acoustic 
propagation.  In  a  layered  atmosphere,  if  the  velocity  is  a  function  of  only  one  space  coordinate 
t  altitude),  the  ray  physics  solution  reduces  to  the  familiar  form  of  Snell’s  law: 


sin  0 
c 


( A  - 1 ) 


where  0  is  the  angle  that  the  ray  makes  with  the  vertical  direction,  cn  and  Q0  represent  a  constant 
reference  velocity  and  ray  angle  at  its  source,  and  p  is  a  constant  parameter  for  each  ray.  Ray  physics 
is  inaccurate  for  quantitative  description  in  such  areas  as  quiet  zones,  ducts,  and  focal  points.  For  the 
predictions  in  this  report,  a  smoothing  process  has  been  employed  that  allows  energy  to  leak  or  diffract 
from  areas  (caustics  and  ducts)  of  high  infrasonic  intensity  to  areas  (quiet  zones)  of  low  intensity 
(see  alternative  expression  for  refraction  anomaly,  A  in  Equation  A-3). 


An  approximation  to  the  propagation  velocity,  c,  in  space  is  obtained  from  wind  and  temperature 
information  vs.  height.  If  there  is  no  wind,  the  sonic  velocity  at  any  height  can  be  calculated  by 

c  -  \Jy kT  =  20.06 v/r  meters/sec,  where: 
c  =  sonic  velocity  T  =  temperature  (/Q 

y  =  ratio  of  specific  heats  k  =  gas  constant  for  dry  air. 

In  the  presence  of  wind, 

c  =  20.06  y/f  +  We,  (A-2) 

where  Wc  =  velocity  of  wind  components  in  meters/sec  in  direction  of  interest.  The  effects  of 
humidity  on  sound  velocity  have  been  considered  negligible  compared  with  other  variables  and 
are  therefore  omitted  in  this  report. 

As  wc  can  sec  from  Equation  A-2,  the  sonic  height  profile  is  obtained  from  a  combination  of 
temperature  and  wind  velocity  gradients,  and  this  sonic  profile  varies  with  direction.  Also,  the  sonic 
velocity  is  anisotropic  since  the  horizontal  wind  components  are  not  the  same  as  the  vertical 
components. 


With  a  description  of  how  the  propagation  velocity  varies  in  space,  wc  can  develop  or  select  an 
appropriate  ray  physics.  Although  the  ray  physics  used  in  this  report  was  developed  by  the  first  author 
from  Equation  A-l  specifically  for  the  BRL  studies,  many  other  investigators  have  developed  ray 
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solutions  for  similar  purposes  (Officer  1958;  Mabry  1962).  The  ray  physics  used  in  this  report  has 
been  adjusted  and  combined  with  semiempirical  equations,  such  as  those  given  in  in  _  -'opendix.  For  a 
static  (or  point)  source,  ray  tracing  is  initiated  at  the  source  for  vanous  elevation  angles  and  azimuths, 
and  these  rays  arc  extended  throughout  the  layered  atmosphere.  This  procedure  is  repeated  for  every 
source  location.  For  the  supersonic  source,  a  directional  ray  tracing  was  developed.  Since  the  ray 
physics  is  not  readily  found  in  the  literature  for  supersonic  sources  following  complicated  trajectories, 
these  techniques  are  briefly  described  in  the  following  paragraph. 


The  sound  field  generated  by  a  source  as  it  moves  supersonically  along  its  trajectory'  can  be 
considered  a  cone  directional  source  that  moves.  At  any  instant,  the  ray  directions  near  the  projectile 
form  a  conical  shape  perpendicular  to  the  shock  cone  as  sketched  in  Figure  A-l,  and  both  the  apex 
angle  (a  function  of  Mach  No.)  and  the  cone  altitude  are  continuously  changing  along  the  trajectory'. 
For  this  ballistic  wave  energy,  the  ray  tracing  part  of  the  program  calculates  and  plots  the  position, 
time,  cusps,  and  intensity  anomaly  which  are  due  to  either  refraction  throughout  the  layered 
atmosphere  or  to  accelerations  of  the  source.  Those  ray  calculations  are  made  only  over  lire 
appropriate  cone;  however,  this  ray  tracing  process  must  be  repeated  all  along  the  trajectory  for  a 
number  of  cones  whose  shape  and  attitude  depend  on  the  speed  and  path  of  the  shell.  Since  the  sonic 
velocities  vary  with  both  altitude  and  direction,  a  three-dimensional  program  is  required.  As  with  the 
point  source,  the  layered  atmosphere  is  determined  from  temperature  and  wind  data  as  described  by 
Equation  A-2.  Also,  as  with  the  point  source,  semiempirical  formulas  must  be  used  along  with  the 
ray  tracing  technique  to  quantitatively  predict  the  intensity  field.  (Of  course,  the  appropriate  formulas 
and  techniques  for  the  moving  source  cannot  be  the  same  as  for  the  point  source.) 


Techniques  predicting  sound  intensity  spectrums  include  the  effect  of  absorption  for  a  range  of 
frequencies. 


/(wO  =  /,(*)  da;“\  (A-3) 

where  l(w)  =  sound  intensity  spectrum  prediction  which  can  be  calculated 
for  a  number  of  frequencies.  If  both  amplitude  and  phase  arc 
required  for  time  histories,  the  complex  notation  for  intensity, 

/,  must  be  used  (see  Equation  A-5).  The  frequency  for  peak 
intensity  will  decrease  with  increasing  distance  from  the  source 
because  of  a,  the  attenuation  coefficient  which  decreases  with 
frequency. 
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Figure  A-l.  Diagram  of  Rays  Generated  by  Projectile  at  a  Point  on  the  Trajectory. 


!j  v.  )  =  sound  intensity  spectrum  known  at  some  point.  For  an  explosive 

source,  IJ w)  is  calculated  by  a  semiempirical  scaling  formula 
at  a  sciectcd  ,  i stance. 

D  =  dive-gen  .  .actor  nom  spherical  or  conical  (cylindrical) 
so'"  icing,  depends  on  the  type  Y  source. 

A  -  i. fraction  anomaly  calculated  from  ra\  density 

.  r/(dr/dQ o)  tanQo\  0<  r/  |  r7  +  ^  r  0l2tanQ  o 
for  a  point  source.  C  is  a  constant. 

r  -  .ait  distance  from  source  to  any  posif  m  in  space. 

y  =  di>.;  ncc  ciong  ray  from  source  to  pose  ion  in  space  for 

prei.  .  ion. 

0Li  =  initial  inci in  toon  '  c-'und  ray. 

a  =  attenuation  coefficient  whicn  varies  with  frequency  and 
altitude,  y.  Since  there  is  some  discrepancy  between 
theory  and  available  empirical  results  in  the  determination 
of  a,  empirical  results  from  the  NASA  rocket  firing 
measurements  have  been  used  with  theory  as  a  guide  to 
develop  a  reasonable  Ijnctional  relationship. 

w  =  2k  F  where  F  is  frequency  in  Hz. 

It  can  be  seen  above  that  a,  the  attenuation  coefficient,  varies  with  the  frequency.  For  the 
frequencies  in  this  report,  the  following  empirical  relationship  was  used  for  dB  attenuation  per 
kilometer  at  the  surface: 

(KdB/km )  =  0.015F  -  0.00001  F2 .  <  -\-4 1 

The  following  equation  is  the  complex  form  of  Equation  A-3.  It  i«  use"*  when  time  histones,  hn, 
as  well  as  spectrums.  I(w),  are  required. 

/(H  )  =  [/  (w)  \F(w)  ! 

L  °  J  L  J  (  A-0  1 

$  \ 

lit)  I.U'i 

where  F(w)  is  the  vector  functional  relationship  between  fa(w)  and  f(w).  f(w)  is  the  vector  power 
spectnjm  at  one  positron  (where  both  amplitude  and  phase  of  /  can  be  plotted  against  frequent  ). 

/"/ Hd  is  the  vector  power  spectrum  at  another  position. 
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The  upward  arrows  in  the  symbolic  system  used  here  (^»)  stand  for  the  Founcr  transforms  ol  the 
signal  intensity-time  histories;  the  downward  arrow  stands  for  the  inverse  Founer  transform  of  signal. 
In  ouicr  words,  by  these  mathematical  transforms  in  Equation  A-5,  indicated  by  the  arrows,  we  can  go 
from  the  signal  time  domain  (as  measured)  to  the  signal  frequency  domain  (as  used  in  the  model)  or 
\  ice  versa.  The  type  of  mathematical  transform  to  be  used  depends  on  whether  the  signal  is 
considered  a  periodic  or  transient  function.  If  the  signal  is  assumed  periodic  (c.g.,  jet  engine),  a  line 
spectrum  is  used  since  the  harmonic  order  assumes  only  discrete  values.  If  the  signal  is  a  transient 
function  (e.g.,  explosive,  initiation,  or  ballistic),  the  Fourier  integral  is  used  (the  period  of  the  series  is 
allowed  to  approach  infinity)  and  the  Fourier  transform  in  this  case  is  represented  by  a  continuous 
spectrum  where  all  frequencies  arc  assumed  possible  instead  of  only  discrete  values.  We  should 
remember  that  the  two  types  of  Fourier  transforms  (line  and  continuous)  can  only  be  compared  for 
relative  values.  In  other  words,  an  upward  adjustment  in  amplitude  for  the  continuous  spectrum  is 
required  if  absolute  comparisons  are  to  be  made  between  the  line  and  continuous  spectrum. 

If  we  refer  to  Equation  A-5,  we  can  see  that  the  vector  functional  relationship  F(w)  can  be 
empirically  determined  from  two  or  more  measurements,  I(t)  and  IJt),  if  the  form  of  this  relationship 
is  known  from  theoretical  considerations.  With  a  good  approximation  of  F(w ),  the  acoustic 
environment  can  be  predicted  at  other  locations. 

F  or  the  supersonic  (moving)  source,  the  BRL  program  uses  an  extension  of  theory 
(Whitham  1950)  to  predict  quantitative  intensities  along  the  rays  in  a  layered  atmosphere.  Two  useful 
results  of  Whitham’s  work  can  be  summarized  in  the  following  two  formulas  (Equations  A-6  and  A-7) 
that  predict  the  excess  pressures  pm  and  time  constant  x  for  the  N  wave  produced  by  a  body  of 
revolution  moving  supersonically. 

Pm  =  K,  Pg  Ps  K2  dlV4  (M2  -  l)1'8  h314,  (A-6 

x  =  Cs‘  1.9  dt1'4  (M2  -  if3  h"4,  (A-7 

where 

pm  =  excess  pressure  of  N  wave, 

K i  =  reflection  coefficient  (usually  close  to  2), 

Pg  =  pressure  at  gage, 

Ps  =  pressure  at  source, 
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K2  =  aerodynamic  shape  factor,  between  0.54  and  0.81, 
d  -  maximum  diameter  of  body, 

/  =  effective  length  of  body, 

M  =  Mach  number, 

h  =  perpendicular  distance  front  gage  to  line  of  body’s  flight  path, 

x  =  time  constant  or  duration  of  /V  wave, 

Cs  =  velocity  of  sound  at  source. 

From  calculations  of  pm  and  x  in  Equations  A-6  and  A-7,  the  time  history'  of  the  ballistic  N  wave 
can  be  determined.  With  the  use  of  the  continuous  Fourier  transform  described  under  Equation  A-5, 
the  power  spectrum  can  be  obtained  from  the  predicted  N  wave  in  the  vicinity  of  the  source,  and  this 
spectrum  can  be  calculated  at  the  greater  distances  with  the  use  of  Equations  A-3  and  A-5. 
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